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SUMMARY 


An analytic model and procedure are described which can be used to 
estimate probability distributions of touchdown conditions, without 
utilizing a Monte Carlo simulation. Example applications are presented, 
including one -step and two-step flare strategies for the space shuttle 
orbiter. The computed probability distributions for the two-step flare 
are compared with the corresponding histograms from an independent Monte 
Carlo simulation of the same situation. The computed distributions are 
very good fits to the Monte -Carlo-generated histograms for pertinent 
longitudinal and lateral variables at three points (including touchdown) 
along the vehicle's final approach trajectory. The model, thus verified 
has more inherent power for assessing the effects of system element 
changes (e.g., vehicle characteristics, flight control laws, ground 
guidance, etc.) than conventional Monte Carlo techniques. 
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INTRODUCTION 


This report describes an approach and landing system model* that can 
be used to estimate distributions of touchdown conditions without utilizing 
a Monte Carlo simulation. The primary purpose of the model is to establish 
a structure containing the system elements, command inputs, disturbances, 
and their interactions in an analytical framework so that the effects of 
changes in the various system elements on control precision and available 
margins of safety can be estimated. The model is applied to the terminal 
phase of flight of the space shuttle orbiter, and example calculations are 
made for the North American Rockwell high cross-range vehicle. These cal- 
culations include a limited application of the model to a one- step, and a 
broader application to a two-step, flare situation. Also included is a 
comparison of the computed two-step results with Bell Aerospace Corporation's 
Monte Carlo simulation results for the same situation. Excellent agreement 
was obtained. 

The report is divided into five sections. Section II contains an 
overview and general description of the analysis technique. Section III - 
presents the one-step flare example application. Section IV contains the 
two-step flare application and a comparison with Bell's Monte Carlo results. 
Section V then follows with conclusions and comments. 


The system. model described herein is intimately associated with particular 
systems analysis techniques and procedures. However, the general name, "system 
model," is used to denote a broader representation of a given situation than 
ty analysis technique." The authors wish to apologize for any 
confusion this semantic choice may createj certainly none was intended. 
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SECTION II 
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OVERVIEW AND GENERAL DESCRIPTION 
OF ANALYSIS TECHNIQUE 

A. OVERVIEW OF ANALYSIS TECHNIQUE 

Rather than using a direct method (such as Monte Carlo) to generate 
touchdown distributions, the final approach phase of flight is broken down 
into two distinct flight segments, and calculations for each segment are 
made separately. The final calculation of touchdown condition distributions 
is then achieved by appropriately combining the separate results. The analy- 
sis is done in this way to make it possible to take advantage of certain 
mathematical simplifications that can be applied only to individual parts of 
the approach calculations. The end result is an alternative to a Monte Carlo 
simulation for estimating distributions of touchdown conditions. 

Basically, the technique is to separate a constant glide-slope phase from the 
flare and touchdown phase. Glide-slope tracking is essentially a stationary 
process involving small perturbations, wherein linearized equations of motion 
are applicable, while the flare and touchdown phase involves larger perturba- 
tions and significant nonlinear effects due to contacting the ground (e • g. 
touchdown only occurs when h is negative and h is zero). 

For the glide- slope tracking phase, linearized equations of motion are 
used to compute rms values of vehicle dispersions due to random gust and 
beam noise disturbances. Because these disturbances can be considered to 
have stationary Gaussian distributions (Ref. 1 ), and because the vehicle 
(plus control system) can be represented as a linear system, the resulting 
vehicle dispersions are stationary and Gaussian. Thus, a set of normally 
distributed (Gaussian) vehicle dispersions can be computed, which represent 
the distributions of initial conditions for the flare phase of flight. 

For the flare phase a nonlinear simulation is run on an analog (or 
digital) computer to obtain a mapping of initial conditions (and the effects 
of wind profiles and shears) onto the ground. The effects of random gusts 
encountered after the start of flare are accounted for by the inclusion of 
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the most significant aspects of a random gust environment near the ground. 
These include a probabilistic ramp wind change during the final few seconds 
prior to touchdown (which will be discussed further later). One more point 
regarding gusts will be mentioned here.' This is that (in our computations) 
the random gusts are effectively "turned off" at the start of flare, which 
is the exact point in the approach that the beam- following part of the flight 
control system is turned off. This coincidence was intentional. If beam- 
following control were allowed to remain on after the (simulated) gust envir- 
onment was discontinued, then the computed gust-induced dispersions would go 
to zero. Clearly, such a procedure should not be adopted. 

A digital computer was used to combine the results of the initial con- 
dition calculations and the mapping of initial conditions onto the ground. 

The computational procedure involves integrations of conditional probabili- 
ties over appropriate variables . (which will be explained in detail later). 

The result is a probability density distribution at touchdown for each 
important variable. To compute the probability of exceeding a given touch- 
down window is then straightforward because areas under tails of probability 
density plots represent exceedance probabilities. 

B. GENERAL DESCRIPTION OF ANALYSIS TECHNIQUE 

^ »• 

Here, it was intended to present a step-by-step list of tasks comprising 
the analysis technique. However, because there are a number of alternative 
groups of steps to be followed at various points in the analysis (depending 
on the particular situation), the step-by-step list gets quite complicated 
if all contingencies are included. To avoid this, we will sacrifice some 
generality in favor of clarity and conciseness. It is unrealistic to imagine 
that we could anticipate all contingencies anyway. Accordingly, we present 
what might be called a "typical" set of steps for the longitudinal situation 
(similar "typical" steps are also required for lateral calculations) as 
follows : 

1 . Simulate complete vehicle, guidance, and control system 
on an analog or digital computer. 

2. Adjust system parameters and initial conditions to obtain 
a nominal trajectory from equilibrium glide through to 
touchdown . 
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3. Compute maximum expected values of dispersions at 
start of final flare. 

© Use linearized, perturbation equations 
of motion for vehicle just prior to final 
flare . 

© Use guidance and control system definitions 
for glide phase. 

© Prom nominal trajectory, determine nominal 
conditions just prior to final flare. 

© Use spectral descriptions of gusts, beam 
. noise, and any other random disturbances. 

© Compute closed-loop transfer functions for 
gust (etc.) inputs. 

© The following steps are suggested to select 
"maximum" gust level. 



0 Ug = °ur x f (altitude) 

a. Pick A = 3.5 a UR = 8 ft/sec [this 
gives P(a UR > A) = 0.0021]. 

b. This gives a gust level (a UR = A) 
to be used to compute cr ^ and a Wg; 
at altitude. 

© Use gust (etc.) power spectra and closed-loop 
vehicle transfer functions to compute power 
spectra for vehicle dispersions, q (q = h e , h € , 
etc.). 

• Integrate power spectra to compute rms values. 

• Pick = 3 C-rj for initial conditions 

at the start of final flare. 

^MAx/ CT q = 3.1 gives P(tj > t^) =" 0.00097] 


one tail = 1/2 P(cx U r > A) 
only 
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© Combine the effects of gusts, beam noise, 
etc., via an RSS calculation to get ranges 
of al 1 initial, conditions at start of final 
flare . 

4. Map initial conditions (at start of final flare) onto the 
ground . 

© Use complete simulation (from l) and apply 
initial conditions one at a time (using 
values up to maximums determined from J>) } 
recording touchdown conditions each time . 

© Apply ramp wind change (representing short- 
time approximation to random gust environ- 
ment during flare) just prior to touchdown. 

Ramp is a horizontal wind change ,, of magnitude 
Aug during time T. With a constant Aug, 
use various durations, T, to determine maximum 
sensitivities for Zpj), X^j), etc. Then record 
touchdown conditions for various Aug magnitudes 
(using times for maximum sensitivities). 

® Apply mean wind profiles of various magnitudes, 
u w , (up to maximum winds) from an altitude of 
1000 ft to touchdown (not just during flare) 
and record touchdown conditions. 

5. Establish touchdown mapping relationships. 

9 Plot touchdown conditions (e.g., X^) vs 

initial conditions, mean wind, and Aug. '*■ 

© Determine best-fit linear equations to plotted 
individual relationships. 

® Combine above equations to get overall equations 
(e.g., Xipp = Xr CDo + AAu w + Hh e + . . .). 

6. Determine important variables. 

© Define "successful" touchdown window. 

© For a worst-case input, determine which 

variables remain well within their respective 
successful ranges. Those that do are most 
likely not important* those that do not are 
the important ones. 

7 . Compute means and rms values (p and a) for the conditional 
distribution of each important variable (Xrp£) is used as an 
example here). 

© Group the terms in the touchdown equations 
as follows : 
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Xtd = - Xqj) 0 + AAu w + Bh e + Cu e + Dh e + EAug 


— ^TDi 


hi 


For a given wind profile, is a constant. 

For a given rras gust level, h-| and EAug are 
independent, zero-mean, Gaussian random 
variables. 


© Use perturbation equations of motion to com- 
pute transfer functions for h-| to Ug, Wg, 
beam noise, and equipment error (etc.) inputs. 

® Use these transfer functions, along with power 
spectra of the disturbances, to compute rms 
values of h-j due to each input. 

© Combine these separate sources of dispersion 
to define a "total" cr^. . 


°hi 

or 



+ 



beam noise 


+ . . . 



( a B /o? )o?. In.? 
' ,v g / a g' xi iwg J 


g 


a?. 

iX ibeam noise 


known from 
gust model 


® Compute the ratio of a Aug to Oug corresponding 
to the best linear fit (using a least square 
error criterion) to a short-time gust history 
as follows: 


° u g 

where: 


cd§T 2 


I 

^24 u^T 5 + 72(2 + cd 0 T ) ( 2 - o^T) - 72(2 + o>oT) 2 e"^° T j 


cdq is the break point of the Ug spectrum 
and T is the duration of Aug determined 
in step 4 


For times less than 10 sec this expression gives 
values which are about 13 percent larger than 
the equivalent values obtained by fitting a 
short-time gust history with a linear fit through 
the initial and end points. The expression obtained 
via this end point technique is given in Ref. 3 as: 

a Au g / a Ug = V4 “ e_m ° T ) 

6 
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Compute conditional mean and rms values. 




a Xpj) = V a hT +e 2 °Au s 




+ 





beam noise 


+ . • • 


8. Integrate conditional distributions over all values of mean 
wind and gust magnitude to obtain "overall" probability- 
density distributions at touchdown for each important 
variable . 

9. Integrate probability density distributions up to the 
"successful" window limits to determine the probability 
of a successful touchdown for each variable. 

10. Combine the individual success probabilities to obtain an 
overall probability of a successful touchdown. 
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SECTION III 


EXAMPLE APPLICATION TO ONE- STEP FLARE 

A. VEHICLE AND TASK DESCRIPTION 

1 . Vehicle Definition 

The example vehicle is the 1 3^C version of the North American Rockwell 
high cross-range orbiter (Ref. 5)- For the example calculations , only the 
longitudinal situation is considered. 

2. Successful Touchdown Window 

Table 1 lists the limiting values of longitudinal variables for a 
successful landing (as used in the example calculations). Note that no 
airspeed requirement is needed because the vehicle is flyable at speeds 
well below the speed that will result in tail scrape. However, excess 
airspeed will affect the stopping distance, X stor> , which is a function of 
the inertial (ground) speed, U. 


TABLE 1 

"SUCCESSFUL" LONGITUDINAL TOUCHDOWN WINDOW 


WINDOW DEFINITION 

COMMENTS 

800 < XqD < 2500 ft 

1550 Ft ±750 ft (insures touchdown 
on runway) 

0 < Z<pp < 6 ft/sec 

3 ft/sec ±3 ft/sec (insures 
acceptable impact at touchdown) 

x stop < 9200 ft 
( X stop = 7000 ft + £%> 

+ 20UJP ft/kt) 

x stop computed stopping dis- 

tance . Upp is perturbation in 
ground speed at touchdown (insures 
stopping on 10,000 ft runway and 
leaves same 800 ft margin as does 

x td) 
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J. Basic Flare Logic 


For the one-step flare an exponential flare path was selected. This 
is accomplished by making commanded sink rate proportional to altitude. 
Details of the selected flare logic will be presented later. 

4. Inputs 

a. Steady Wind 


The steady" wind profile of Ref. 2 was used in the example cal- 
culations. This results in a profile whose magnitude and direction 
are determined by random selections from given distributions. Thus, 
for any given approach and landing, the profile is fixed, but from 
one approach to the next the profile changes. A sample profile is 
shown in Fig. 1. To obtain any other profile it is only necessary 
to scale up (or down) the wind magnitude. Conveniently, any parti- 
cular profile can be completely determined by specifying the magnitude 
at a given reference altitude. For our purposes a wind reference 
altitude of 10 ft. was selected. At. this altitude the wind magnitude 
varies from a 10 kt tail wind to a 25 kt head wind. 

By combining the distributions of magnitude and direction given * 
in Ref. 2, a distribution of the wind component parallel to the runway 
can be determined. A rough plot of this distribution is shown in 
Fig. 2. It is presented to give a feel for relative likelihoods of 
the head and tail winds that are used. From the plot it is seen that 
the most likely longitudinal wind condition is a head wind of about 
8 kt. As a consequence, a "nominal" flight condition will be defined 
that includes an 8 kt head wind (measured at an altitude of 10 ft). 

b. Random Gusts 

In addition to the steady wind profile, a turbulent wind environ- 
ment is also considered. The model for this turbulent wind is given 
in Ref. 2. Basically, the model provides fop random gusts whose rms 
level has a Rayleigh distribution (with a characteristic speed of 
2.5 ft/sec), and whose frequency content is a function of altitude. 
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Wind (kts) 


Figure 1 . Wind Profile Associated with 
10 Kt Wind at 1 0 Ft Altitude 
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Figure 2. Approximate Distribution of^Head and Tailwinds (Measured at an Altitude of 10 Ft) 


(Table 2 gives the gust spectral shapes for an altitude of 1200 ft . ) 
Further, the model provides for the ratio of the rms level of ver- 
tical gusts to the rms level of longitudinal gusts to also be a 
function of altitude (although the individual gust components are 
not correlated). At 1200 ft the ratio is 0.88. 

TABLE 2 

GUST POWER SPECTRA FOR AW ALTITUDE OF 1 200 FT* 


u g 

Ku 

$U g U g o 2 + 0.11 

w g 

k2(u>2 + 0.66) 

w g w g (o>2 + 0 . 1 8)2 


* Values of K u , K w are selected 
from a Rayleigh distribution 


For our calculations, these random gusts are applied to the 
vehicle only prior to the start of flare. However, because a random 
gust environment exists all the way to touchdown it is necessary to 
account for any significant gust effects during flare by adding back 
the most important gust effect, which is a ramp change in horizontal 
wind during the last few seconds prior to touchdown. The rest of the 
random gust environment (during flare) has very little effect on the 
touchdown conditions. This gust-produced ramp change in wind is 
actually a random wind shear, as discussed in the next paragraph. 

c. Random Wind Shear 


As mentioned above, the random wind shear is the most significant 
aspect of the random gusts during flare. Although its primary effect 
is to produce hard landings, the effects on other touchdown variables 
were not neglected. A hard landing generally arises because an 
unanticipated large decrease in vehicle airspeed (due to a tailwind 
gust or a decrease in headwind) occurs when the vehicle is close to 
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the ground. The decrement in airspeed produces a decrement in 
lift (even though the angle of attack may be increasing ), which 
in turn produces an increased sink rate, which then results in an 
almost immediate touchdown. If the- random shear is of the opposite 
sign (i.e., an increase in airspeed) then the flight path is shal- 
lowed out somewhat, but nothing "interesting" happens. 

The random shear is simulated as a ramp change in horizontal 
wind (superimposed on the shear due to the steady- wind profile) 
starting at an. altitude of approximately 10 ft. This results in 
about 3 sec of shear duration prior to touchdown, which is consis- 
tent with the results of Ref. 3 (for conditions leading to hard 
landings). The magnitude of this wind change is given in Ref. 3 
as a zero-mean, Gaussian random variable with an rms value defined 
as 

CT Aug - % V 2 0 - e~ r/La ) 0) 

where r is the distance flown (during the Aug), and Ly is the scale 
length for longitudinal gusts. However, this representation is one 
in which a linear fit merely connects the initial and final points 
of a sample time history of Ug. As noted above, in the listed 
analysis "steps, " a better representation of a short-time gust 
history is obtained by fitting the entire sample, rather than just 
connecting the end points. The corresponding expression for o Al1 

O 

using a least squares fit is: 


g Au ct = -TT% Js^T 5 + 72(2 + c^T)(2 - cd d T) - 72(2 + o^S) 2 e ~™° T (2) 
s cdoT^ ’ 

where cuq is the breakpoint of the Ug spectrum and T is the duration 
of the sample. (Note that r/I^ used earlier is identical to cd q T used 
here.) This preferred representation results in c f/w being about 
13 percent larger than the value obtained for the. end point fit (for 
times up to 10 sec).. 
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d. Other Error Sources 


Glide slope beam noise , measurement errors, pilot remnant 
(for manual control), and other sources of errors applicable 
in general, have not been included in these first simple 
illustrative example calculations. However, they are included, 
as appropriate, in the later two-step flare calculations. 

5. Preflare Plight Condition 

The steady-state flight condition prior to flare is defined in Table 3- 


TABLE 3 

PREFLARE FLIGHT CONDITION 


Flight path angle 

—10 deg 

Airspeed 

302 kt 

Angle of attack 

5.3 deg 

Nominal wind at 
800 ft altitude 

- 16.5 kt 

(headwind) 


6. Preflare Plight Control System 

Figure 3 shows the inner- and outer-loop control system block diagrams 
for preflare glide slope tracking. Table 4 contains a list of the dimen- 
sional stability derivatives that are appropriate for the preflare steady- 
state flight condition. The associated closed-loop Bode plots are presented 
in Figs. 4 and 5 . 
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a) Inner Loop SAS Gives Augmented Vehicle 





b) Outer Loop Control System 

Figure 3- Closed- Loop Block Diagrams for Preflare Flight Control System 
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TABLE 4 


DIMENSIONAL DERIVATIVES FOR PREFLARE FLIGHT CONDITION 


x u 

= — 0 . 0206 1 / sec 

z & 

°e 

= - 242 . ft/sec 2 /rad 

*w 

= —O.OO654 l/sec 


= 0 

x & 

°e 

= — 40.7 ft/sec 2 /rad 

“w 

= -O.OO579 1 /ft- sec 

Z-u 

= -0.126 1 /sec 

Mq 

- -O.579 l/sec 

Zw 

= —1 .28 1 /sec 

M6 e 

= ~ 3 . 4 l l/sec 2 


7 . Calculation of RMS Dispersions Prior to Flare 


The Table 4 values, and the control/guidance logic represented by the 
block diagrams of Fig. 3, were used in standard (e.g., Ref. 4) linearized, 
closed-loop perturbation equations of motion to obtain transfer functions 
for gust inputs. (The effect of the nominal wind condition was accounted 
for in the equations of motion by utilizing airspeed in the terms arising 
from aerodynamic forces, and ground speed in the terms arising from inertial 
accelerations. Although this is not a large effect, it was done to minimize 
the introduction of errors in the calculations.) 


Using the gust spectra (given earlier) with the gust input transfer 
functions (given in the Appendix) leads to power spectral density expres- 
sions for each variable. Because the area under a power spectral density 
plot is the variance for the input considered, and because variances 
due to independent Gaussian inputs can be added directly, the overall 
rms value for any variable can easily be found by taking the square 
root of the sum of the variances from all of the random inputs. For our 
situation this becomes (using q as a generic example variable): 



JL 

Jo 

u g 


Vs 4 " + 


f 

Jo 


JL| 

Wg 


^WgWg^ 


1/2 


(3) 
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Expressions such as this were used to compute rms values for vehicle 
dispersions just prior to flare. An rms horizontal gust level of 2.04 ft/sec 
(which is the most probable gust level according to the given Rayleigh dis- 
tribution at 1200 ft altitude) was used to compute rms dispersions in speed, 
sink rate, and deviation from the glide slope beam. The results are: 

cr u = O.837 ft/sec 
cf£ = 1.25 ft/sec 

= 2.47 ft 

Further discussion of rms dispersions is deferred until the equations for 
touchdown conditions have been presented. 

8. Flare Control System 

The flare control system used in the example calculations is a constant 
gain system that is switched in at the flare initiation height of 800 ft. 

It was possible to use a fixed-gain system because the airspeed at touch- 
down is sufficiently high that the 6 c /h e ratio remains acceptable. (However, 
an airspeed loop was used to effect a continuous trim change as speed bleeds 
off.) The outer-loop block diagram of the overall control system (including 
the switching required to go from the preflare system to the flare system) 
is given in Fig. 6a. As indicated in Fig. 6b, the landing gear is lowered 
as a function of sink rate at 4000 ft altitude to control the landing speed 
(in the presence of winds) without making large changes in the flare trajectory. 

Feeding forward an open-loop pitch or pitch rate command (an early variant 
of the flare system) to minimize deviations from an exponential trajectory did 
not significantly affect touchdown performance. What was really needed was 
an additional pitch-up command just before touchdown to compensate for the 
loss of lift due to speed bleeding off. But, rather than doing this with an 
open-loop feedforward command, it was done via an airspeed loop closure (during 
flare only). Although this loop is unstable in the classical sense (because 
a decrease in airspeed results in an increase in pitch attitude which causes 
a further decrease in airspeed, etc.), in this situation touchdown occurs 
before the divergence has a change to build up. The airspeed loop has two 
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Figure 6b. Gear-Down Altitude as a Function 
of Sink Rate at 4000 Ft Altitude 



distinct advantages. It provides the desired increase in pitch rate near 
touchdown (when airspeed bleedoff cancels the expected lift increment from 
any angle -of -attack increase demanded by the sink rate error control system), 
and it provides windproofing against shears that otherwise might produce 
hard landings. That is, if a sudden decrease in headwind occurs near touch- 
down, the sink rate increase (due to the airspeed decrease) is averted by 
pitching up to increase the angle of attack (and maintain lift) before a 
sink rate error can develop. 

The magnitude of the gain in the airspeed feedback was determined as 
follows. The lift equation is 

L = 

With lift equal to weight, 

a 




2W 

pu|sc L , 


‘a 


(b) 


( 5 ) 


Then, 


da 

dU a 


— 4w 


PUI SC L, 


'a 



Therefore, 


Aa 1 


—£’ a o 

U a 


AUp 


For relatively short times. 


giving 


A9 = Aa 


— 2ao 

A8 = AU„ 

U a a 


( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 
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( 10 ) 

( 11 ) 


( 12 ) 

This value of gain gave good results. A higher value was tried (to compensate 
for lags in pitch response) but it resulted in worse sink rate control. 

B. ANALOG SIMULATION TO MAP INITIAL 
CONDITIONS ONTO THE GROUND 

1 . Equations Used 

The equations mechanized on an analog computer were the "total" equations 
of motion, as opposed to only perturbation equations. Further, the dynamic 
pressure was one of the variables (rather than just using an average value). 
Based on the specifically applicable Ref. 5 data, the nondiniensional sta- 
bility derivatives were held constant during the flare. The equations that 
were modeled and the coordinate system definitions are presented in Fig. 7 . 

In this figure the subscripts I, a, w, and B refer to inertial, aerodynamic, 
wind, and body, respectively. The rest of the symbols are standard, and are 
not redefined here. Figure 8 presents an analog mechanization diagram to 
show how each variable was computed. 

2. Results of Mapping Initial 
Conditions (IC’s) and Wind Effects 

Figures 9 a through 9 * 1 present the results of mapping initial conditions 
at the start of flare onto the ground, one at a time. Also included are 
the effects of wind profiles of various magnitudes (u w is the wind speed at 
an altitude of 10 ft) and wind shears near the ground (Aug). On each plot 
is a linear fit to the data, and the mathematical expression for the linear 
fit. It can be seen that linear fits are very good approximations to the 
data, with the possible exception of the fit for negative values of 


For our flare situation. 


cto = 11 deg 

= 330 ■ ft/sec 


Ua = 


Therefore, 


A8 

AU. 


a 


-o-o* wh 
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C Li = C I 0 + CLa^B 
C L “ ^Li + ^Lgg^e 
CD = Cd q + ^Cl? + CD5 e Se + Cife g 5g 

x 

Cm = C roo + C^a-g + Cmg^&e + ^y a C m qq 


U I - u a cos °w + % cos 7l = u a + % 

0 = 7 a + «B 

71 = 7 a + Qw 


/- u„ sin 71 \ 

. -1 /” % sin 7a\ 

. “ Uwla 

\ u a j 

~ S1 M uj ) 

% 


~ Ua 

X = mV i = QS £- Cd cos + C L sin c^] - mg sin 7 j 

Z = — mUj7l = Q s t” C L cos “w - C D sin + m § cos 7l 
M = Iy0 : QSc’Cjj! 


Figure 7- Equations Modeled and Coordinate 
System Definitions 
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Figure 8. Analog Mechanization Diagram 
























Figure 9- Mapping of Initial Conditions onto the Ground 
a) due to h £ , u £ 
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Figure 9. Mapping of Initial Conditions onto the Ground (continued) 

b) AX^ due to Au g , Au w 
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Figure 9* Mapping of Initial Conditions onto the Ground (continued) 

c) AZyp due to h e , u £ * 
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Figure 9. Mapping of Initial Conditions onto the Ground (continued) 

e) £LJ td due to h £ , u £ 
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Figure 9. Mapping of Initial Conditions onto the Ground (continued) 

f) AUjj) due to Aug, Au-w 
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Figure 9- Mapping of Initial Conditions onto the Ground (continued) 

g) A0rjj) due to h e , u e 
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Figure 9 . Mapping of Initial Conditions onto the Ground (concluded) 

h) ASrp-Q due to Au^, Au w 
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Aug. However, this side of the curve is not of great significance because 
it doesn't tend to produce hard landings. Superposition of the effects of 
several inputs was verified as being valid. A summary of these linear fits 
is presented in Table 5 (which gives equations for the touchdown conditions 
as functions of the initial conditions and winds). Also included is the 
expression for stopping distance (X s t 0 p)* For the system considered, beam 
deviation (height) errors at the start of flare (h e ) have an effect only on 
and X sto p. The effect is a translation in X by an amount equal to 
h e cot 7 (because there is no height error at the start of flare — only an 
X error). 


TABLE 5 

TOUCHDOWN DISPERSION EQUATIONS 


X TD 

Zw 

UTD 

0TD 

x stop 


= 1550 + 5-°Au w + 5-67h e + 35*3u e + 77*5h e - i5.5Au g 

= 2.83 - 0 . 03 Auw - O.Ol54u e + 0.0475h e + 0.0987Au g 

= 297.5 + 2.89Au w + 0.193 u € - I.l9he + 0.275Au g 

= 11.86 - 0.052Au w - 0.01 l8u e + 0.076h e + 0.0445Au g 

= 7000 + AXtd + 33.25*AUpd 


-J 

0 

0 

0 

+ 

1 01 Au w + 5.67h e 

+ 4l ,7u e + 

37-9h e - 

Units: 

Au w , 

kt 

^TDj 

ft 


u e . 

ft/sec 

• 

Z TD, 

ft/sec 



ft/sec 

U TD, 

ft /sec 



ft 

e TD> 

deg 



ft/sec 

x stop> 


*This number is based on the assumption of a constant deceleration of 
10 ft/sec 2 beginning 3.5 sec after touchdown. . 
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3* Determination of Important Touchdown Variables 


Now that the equations for the touchdown conditions have been determined, 
we need to find out which, touchdown variables are important and which (if any) 
may be ignored. This is done by checking the variables (if any) in Table 5 
that remain well within their respective acceptable regions when a worst-case 
input is applied. The worst-case input is approximated by the maximum wind 
magnitude and a 3 ° gust level. RMS values for the touchdown variables are 
found via a technique described in the next subsection. The results are: 


° X TD 

= 324 ft 

°Ztd 

= O.85 f t/se c 

ct 0 TD 

= 0.47 deg 

y 

stop 

= 183 ft 


In Table 1 the limiting acceptable values for each variable were 
given. Combining these with the nominal values from Table 5 gives a margin 
for each variable. These are 


Mx td - 
M ZpD = 
M 0 T d = 

Mx stop = 


'750 ft 
3.2 ft/sec 
3.1 deg 
2200 ft 


Decreasing these margins to account for the effect of £u Tr = ±18 kt 

W' y 

and then dividing the modified margins by the rms values leads to a rough 
idea of the likelihood of each variable exceeding its margin (for the 
given worst-case input). Thus, 


a *TD 


$60 

32 ^ ~ 


2.0 
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3.2 


t 


MZTD 

. 2.7 . 

CT ZtD 

0.85 

Mgrjrp 

i 2.2 

ff 0TD 

0.47 

^stop 

. 382 . 

°^stop 

183 


For a Gaussian distribution (which applies to this case because it has a 
given u w and a 15 „) the probability of exceeding the margin when M 'fa = 2 

’ 1 ' [ O 

is about 0.046. For M'/cr = 3 the probability drops to about 0.0027; an d- 
for M'/a = 4 it drops to about 0.00006. Thus the probability of the pitch 
attitude at touchdown exceeding its maximum acceptable value is insig- 
nificant compared with the probability of Xjp being outside the ±750 ft 
acceptable window. 

Although all of the probabilities may be small, here we are only con- 
cerned with relative likelihoods. We are merely determining which variables 
to consider further. Based on the above results, the conclusion is to 
consider Xpp, Zpp, and X s . 

C. COMPUTATION PROCEDURE 

The computations are best explained by considering an example case. 

For the example chosen (X^), the dispersion equation given in Table 5 is 
of the form, 

XjjTQ = + AAu w + Bh e + Cu g + Dh^ + EAUg 03) 

where A, B, C, D, E are known constants. The right side of this equation 
can be broken into three parts using the following definitions. 



%>! 

= ^TDq + AAu-w 

0>0 

and. 

hi 

= Bh e + Cu e + Dh £ 

(15) 
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Xijdi + lii + EAu.g 


06 ) 


Thus, 


X TD = 


The logic behind grouping the terms this way is simple. For a given wind 
profile, XpD] is a constant. For a given Gaussian gust environment, both 
h^ and EAug are Gaussian (with zero means) and independent. Further, both 
h-j and EAug have rms values* that are proportional to Og. t Therefore, for 
a given wind profile and rms gust level, X^ has a Gaussian distribution 
with a mean given by 


^XpD x TD-j 


(IT) 


and an rms value given by 


°Xtd 


or, 




TD 


\/ °h^ + E?0 Iu g 
B 2 a g 


08 ) 

09 ) 


With the wind profile and gust level specified, the distribution of 
Xijijj is actually a conditional distribution. As noted above, this condi- 
tional distribution is Gaussian. However, as will be shown later, the 
overall distribution of Xqtq is not Gaussian, This is pointed out here 
to emphasize the fact that Gaussian gusts and a set of linearized touchdown 
relations do not necessarily produce a Gaussian distribution for touch- 
down dispersions. One reason for this is that the rms gust level is not 
a constant, but is itself a random variable (with a Rayleigh distribution). 

With a Gaussian conditional probability density distribution for Xpp, 
it is a simple matter to compute the conditional probability of a success- 
ful touchdown. The unshaded area in Fig. 10 represents this conditional 


*Because hi is a linear sum of perturbation variables, an auxiliary 
equation can be used in the matrix of perturbation equations of motion to 
obtain transfer functions for hi to ug and wg inputs directly. Then ohi 
can be found (for a given cr U g) via the technique described earlier. 

*In this derivation cr U g will be called Og to simplify the notation. 
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P(XtD I U w , Og) 



Figure 10. Conditional Probability Density Distribution for Xtd 

probability. Thus the conditional probability of a successful touchdown 
(given u w and a g ) is:* 

P(X TD € X S |u w , cr g ) = 1 - F Xl ~ Fx 2 (20) 


This conditional probability can be multiplied by the probability 
density distribution of u w , and integrated over all u w to give another 
conditional probability; one that depends only on a g . (A rough plot of 
the probability density distribution of u w was given in Fig. 2.) Thus 
the conditional probability of a successful touchdown (given cr g ) is. 


P(X T d 


/ OO 

P(*TD 


With the conditional probability a function of only one variable (a ) 

O 

it can be plotted to show the effect gust level has on touchdown conditions. 
Such a plot looks like that shown in Fig . 11. 


*The notation used in this equation is read as follows. "The Proba- 
bility of Xpp belonging to the successful range of Xpp, given u w and o g , is 
1 — Fx-| — Fx 2 «" The symbol Xg is used to denote the successful range of X^. 


TR- 1007-1 


38 




Figure 11. Conditional Probability of a Successful Touchdown 
as a Function of RMS Gust Level 

It remains now only to multiply this conditional probability by the 
probability density distribution of cr g , and integrate over all a g to obtain 
a number for the overall probability of a successful touchdown (i.e., be 
within X-pp limits). This is shown in the next equation. 


P(X TD e Xg) = J P(Xtd « Xgjcr g )p(ffg)dag (22) 

Probability of , 
a successful 
touchdown 

The above explanation considers only one way to carry out the required 
calculations. It is also possible to carry out the various integrations in 
a different order to end up with the last integration over X. Then the 
overall probability density distribution of X T p becomes available (rather 
than a conditional version) as shown in Fig. 12, and discussed below. 

D. TOUCHDOWN DISTRIBUTIONS 

Figure 12 shows a plot of the overall probability density distribution 
for Xtd* The cumulative area under the curve (integrating from right to 
left) represents the probability of exceeding any given value of X W . This 
exceedance probability is shown plotted on probability paper in Fig. 13 to 
show that the distribution of Xpp is not Gaussian. Gaussian distributions 
plot as straight lines on probability paper. • In particular, it can be seen 
in Fig. 13 that a Gaussian distribution would underestimate the probability 
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P (X TD >X TD ) Percent 



Xt D (ft ) 


Figure 13* Exceedance Probability of X^d 

4l 
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P (X TD < X TD ) Percent 




of very large dispersions from the mean. This is an encouraging result 
because measured data commonly give more large dispersions from the mean 
than would be predicted by a Gaussian distribution that matches the mean 
and variance . 

A plot of the probability density distribution for sink rate at touch- 
down is presented in Fig. l 4 , and the associated exceedance probability is 
in Fig. 15. Plots of stopping distance (X s ^ 0 p) probabilities are not 
included because the probabilities turned out to be significantly smaller 
than those for Xpj) and Zpp. This was a surprise since the value of U'/o 
was not large. It turns out that the computed value of M*/a is misleading 
because the probability of the 18 kt Au w that was used is nil. Had a 
0.3 percent value of Au w been used, then M'/cr would have been about 4 . 3 > 
which is a more accurate representation of the relative importance of Xg-fcop* 
In any event, this points out the fact that X s top is sensitive to tailwinds 
for the example system. 

E. PROBABILITIES OF EXCEEDING TOUCHDOWN WINDOWS 

From the distributions of touchdown conditions it is found that the 
acceptable touchdown windows are large compared to the expected touchdown 
dispersions. Clearly, this is desirable for an acceptable system. From 
Fig. 13 the probability of exceeding the allowable Xpp window is found to be: 

pCxrjD < 800 ft) + PCxpp > 2300 ft) 1 600 X 10 -6 

The probability of exceeding the touchdown sink rate limit of 6 ft/sec does 
not fall within the bounds of the plot in Fig. 15. However, the probability 
was computed to be: 


PCZpp > 6 ft/sec) = 3 x 10" 6 

These numbers make it obvious that short and long landings are much more 
likely than a hard landing for the example system. Therefore, a system 
improvement could be made by reducing the Xpp exceedances at the expense 
of larger Zpp dispersions. 
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Figure 14. Probability Density Distribution for Sink Rate at Touchdown 
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Z TD (ft/sec) 


Figure 15* 


Exceedance Probability for Z^p 
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A few comments on the primary causes of the touchdown dispersions are 
pertinent here. Using the equations in Table 5 and the computed vehicle 
dispersions at the start of flare indicates that AZpp comes mainly from 
Aug, while AKpp depends strongly on sink rate at the start of flare, and 
does not depend very strongly on Aug. Thus, in addition to defining the 
main sources of dispersion, this shows that the important touchdown dis- 
persions will be largely uncorrelated. 

F. COMMENTS AND CONCLUSIONS REGARDING 
ONE -STEP FLARE EXAMPLE 

® Although not mentioned earlier, it was found that the 
speed during preflare glide must be considerably higher 
than the speed for maximum L/D in order to provide an 
adequate margin for increasing L/D when a large head- 
wind is present. In terms of flight path angle, our 
preliminary calcinations indicate a descent angle of at 
least 9 beg is required. 

® To achieve reasonable touchdown conditions in the presence 
of winds, some form of speed control is required. Varying 
the height for lowering the landing gear as a function of 
expected speed error at touchdown (similar to what was 
done for the example system) is a simple way to reduce 
speed errors to acceptable values. 

© The variables that had the most effect on touchdown con- 
ditions for our example system were head- and tailwinds, 
and speed and sink rate prior to flare. 

© The computed distributions of touchdown conditions showed 
the tails to be significantly larger than those of equi- 
valent Gaussian distributions. As a consequence, equiva- 
lent Gaussian distributions would appreciably underestimate 
the probability of large disperions from the mean. 

• For large tailwinds, the probability of running off the 
far end of the runway goes up very rapidly. The reason 
the computed probability of exceeding the acceptable 
stopping distance was not large is that the given wind 
model provides a negligible probability for tailwinds 
over about 8 kt. 

© By simulating the entire nonlinear situation on the analog 
computer, and then using the empirical (measured) results 
in the way discussed, it was possible to make use of very 
powerful linear analysis techniques that otherwise would 
not have been applicable. 



SECTION IV 


TWO-STEP FLARE APPLICATION 


This section presents the results of an analysis of Bell Aerospace 
Corporation's two-step flare and decrab control system for the North 
American Rockwell high cross-range space shuttle orbit er vehicle. Figure 1 6 
presents a side view of the two-step flare trajectory, showing the initial 
flare (from an equilibrium glide), the constant flight-path angle glide 
phase, and the final flare t-o touchdown. With a sufficiently long glide 
phase (of the order of 15 sec or more) it is not necessary to consider the 
effects of the initial flare on the distributions at the start of final 
flare (because all transients have had time to die out) . This is the case 
with Bell's guidance strategy. Thus we need only consider the glide and 
final flare flight segments in our analysis (except when determining the 
effects due to various mean wind profiles) . This reduces the two-step flare 
analysis to an effective one-step flare. Therefore, rather than repeating 
a description of all the steps in the analysis, only the results will be 
presented. 



Figure 16. Side View of Two-Step Flare Trajectory 
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Included as Fig . 1 7 are the longitudinal and lateral block diagrams for 
Bell's control system (Ref. 6). Figure 18 then presents the results of 
mapping the conditions just prior to final flare onto the ground. Discussion 
of Figs. 17 and 18 is limited to only the following two comments because 
the results are the same as were found previously for the one-step flare. 

First, it is pointed out that linear fits to the touchdown conditions are 
very good representations of system sensitivities . And second, the variable 
Ahpp. is included to show the effect of perturbations in the altitude at 
which final flare is initiated. One further comment should be made here. 

This is that Bell ran their Monte Carlo runs in two separate series; one 
with a scanning beam landing guidance system, and the other with a tracking 
radar for landing guidance information. In most variables both systems 
gave essentially the same results. For these variables we took the liberty 
of combining the separate results to obtain a single combined result. However, 
for those variables that showed a significant difference for the two systems, 
we accordingly made separate comparisons . 

The primary results to be discussed in this section are the comparisons 
of Bell's and STI's distributions of pertinent variables just before final 
flare, just before decrab, and at touchdown. It was decided to apply the 
approach and landing model at these discrete points so that any significant 
differences in touchdown distributions ( if they occurred) might be traced 
back up the approach trajectory to determine their^ origin and causes. How- 
ever, our preliminary calculations showed no great differences with Bell's 
results. Thus, the primary purpose of the multiple -point comparisons had 
been obviated. As a consequence, we expended only a limited effort on the 
computations of the distributions prior to touchdown, which resulted in 
some of the variables being compared via plots of the distributions, and 
others being compared on the basis of the rms values of the dispersions . 

A. COMPARISONS JUST PRIOR TO FINAL FLARE 

^FF 


On Bell's Monte Carlo runs a one-second sampling time was used for 
printing out the pertinent variables . Because the ground speed just prior 
to final flare is about 325 ft/sec, a one-second sampling time introduces a 
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Figure 17b. Bell Shallow-Glide Guidance Model 
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Figure 1 Yd. . Eell Lateral Control System (Prior to Decrab) 
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Figure 18. Mapping of Initial Conditions onto the Ground (continued) 

b) AXipp due to Aug, Ahpp 
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Figure l8. Mapping of Initial Conditions onto the Ground (continued) 

c) AZpp due to h 6 , u £ 
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dispersion in the printed-out values of X FF that can be up to about ±1 62 ft 
( = ±l/2 X 325). For a uniform distribution of sampling times (about the 
nominal) , an rms Xpppop of about 95 ft is thereby introduced into the dis- 
tribution supplied to us . (Actual X-errors are produced by equipment 
imperfections , random gusts and mean wind variations, and by beam bias 
and noise) . 

In spite of this situation, we can obtain a rough estimate of the size 
of the measured X FF dispersions by using a root -sum -square (RSS) calculation 
to combine the rms sampling errors with our estimate of the actual rms Xpp 
dispersions. This calculation gave us rms values of l8l ft and 1 67 ft, 
respectively, for the scanning beam and tracking radar systems, which 
compare very favorably with the rms values of 1 9 6 ft and 1 83 ft supplied 
to us from Bell's data. 

Upp 

The distribution of inertial speeds just before final flare is due almost 
entirely to the variation in the mean wind. Thus, we can easily make an 
estimate of the expected distribution. According to Bell's prediction 
scheme (indicated in Fig. 17b and detailed in Ref. 6), the ground speed 
at the start of final flare is modified by an amount equal to half of the 
mean wind speed. Therefore, our estimate of the effect of the mean wind 
on Upp is as shown in Fig. 19. Also shown in Fig. 1 9 is a combined histo- 
gram of Bell's measured distributions. As seen in the figure, the comparison 
is quite good . 

Zpp 


Here is another case wherein a significant contribution is due to mean 

winds . Our computed dispersion due to random gusts has an rms value of about 

0.6 ft/sec. The effect of mean wind variations has been estimated at 

0.86 ft/sec rms. (This estimate is based on the assumption that the ratio 

of to My , due to mean winds, is the same as the ratio of a,. to 

npp ^FF' u w?F 

Mupp • ) The combined RSS value for these two dispersion sources is 1 .05 ft/sec, 
which compares very favorably with the 1 .15 ft/sec average for Bell’s two 
systems . 
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Estimate of probability 
density of U FF due to 
mean wind at 130ft altitude 


Bell's histogram 
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Figure 19. Comparison of Bell Distribution of Inertial Speed 
at Final Flare with Estimated Inertial Speed Due to Mean Wind Only 
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■^FF 


A comparison of rms levels of ypp can also be made by RSS'ing our 
separate results . For this variable the primary contribution is due to 
beam bias errors. Thus, separate comparisons for each guidance system are 
appropriate here. When the RSS calculations are made, we get 7-5 ft for 
an estimate of the rms ypp for a scanning beam, and 4 .8 ft for a tracking 
radar. Bell's corresponding values are 7-7 ft and 4.0 ft, respectively. 
Thus, a very good comparison is again obtained. 

It is noted that the effects of mean winds on the lateral position (and 
rate) distribution( s) just prior to final flare are much less significant 
than they are for the longitudinal variables discussed earlier. 

ypp 

Our calculations show that the distribution of lateral position rate 
is relatively insignificantly affected by mean winds and equipment errors . 
Thus, we expect the major contribution to ypp to be from random gusts. 
Figure 20 shows a comparison of our computed distribution of y F p (due to 
gusts alone) with Bell's histogram of y F p (due to all disturbance sources). 
In the figure it can be seen that the comparison is excellent, indicating 
that the major source of ypp dispersions is indeed from random gusts . 

CPFF 


A comparison of Bell's and our computed bank angle distributions just 
before final flare was also obtained, and is shown in Fig. 21 . It is seen 
that the comparison is again exceptionally good. 

B. COMPARISONS AT DECRAB 

The comparisons of lateral variables just prior to decrab are almost 
identical to those presented above. Therefore, we will not repeat the pre- 
sentation, but rather we will go directly to the touchdown comparisons . 
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Bell's Number of Occurrences per.5ft/sec Interval 


p(y FF ) 



y FF (ft/sec) 

Figure 20. Comparison of Bell Distribution of y at Final Flare 
with Computed y^p Due to Random Gusts 
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Figure 21. Comparison of Bell Distribution of cp at Final Flare 
with Computed, cppp Due to Random Gusts 
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C. COMPARISONS AT TOUCHDOWN 


^TD 

Because the scanning beam and tracking radar systems give rms touchdown 
dispersions that are noticeably different, we have made separate comparisons 
in Fig. 22, as well as a combined comparison in Fig. 23. Excellent corre- 
spondence is evident (separately and combined) . Figure 24 presents a 
comparison of cumulative distributions of Xpp which indicates that the 
calculated distribution is the sane as Bell's, except that it is shifted 
about 70 ft down the runway. The reason for the shift is not known at this 
time. However, it could be caused by any number of things, such as a 
difference between c .g . and glide slope antenna locations, or by an equip- 
ment bias of about 8 in. (corresponding to the change in c .g . height due to 
landing gear compression) . 



The comparison of touchdown speed distributions is shown in Fig. 25- 
Again, excellent agreement is found between Bell's Monte Carlo results and 
the analytical-empirical results. As an item of interest, it is noted 
that the double peaking in the computed distribution is believed to be a 
result of the particular model used for the mean wind variations . Figure 26 
shows that Bell's histogram for the scanning beam guidance system also 
exhibits such a characteristic . 

Z TD 

The distribution of sink rates at touchdown is the only area in which 
any difference was found between Bell's and the calculated results. Although 
this difference is small. Fig. 27 does show that Bell obtained a slightly 
wider distribution than we computed for sink rate at touchdown . All in all, 
the fit is actually pretty good, but we mention the difference because it is 
the only one we encountered in all of our comparisons . 
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Bell's Number of Occurrences Bell's Number of Occurrences' 

per 250ft Interval per 250ft Interval 



P (X-fp) 




p(XjD ) 
_l_ 
ft 


Figure 22. Comparison of Bell’s Monte Carlo Generated 
Histograms of Xpp with Computed Distributions 
(for Two Guidance Systems) 
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Figure 23. Comparison of Bell's Monte Carlo Generated 
Histogram of Xpp with Computed Distribution 
(with Data for Both Guidance Systems Combined) 
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Figure 25. Comparison of Bell's Monte Carlo Generated Histogram of -with the Computed 
Distribution (with Data for Both Guidance Systems Combined) 
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Figure 26. Comparison of Bell's Monte Carlo Generated Histogram of with the 
Computed Distribution (Using Data for Scanning Beam Guidance System) 




Figure 27. Comparison of Bell's Monte Carlo Generated Histogram of Z TD with the 
Computed Distribution (with Data for Both Guidance Systems Combined) 



^TD TtD 


For large crosswinds, Bell’s decrab system causes the vehicle to finish 
its alignment maneuver several seconds prior to touchdown. This leads to 
rather large drift rates and lateral offsets at touchdown, as evidenced in 
Bell’s touchdown distributions. However, our original model was not designed 
to include such long decrab times. Therefore, a modification to the model 
would be required to cover this unanticipated situation. But, considering 
the tradeoffs involved in making the required changes, it was decided not 
to expend the effort to modify our calculation technique to cover this 
situation. As a result, we did not compute touchdown distributions of y and 

o 

y to compare with Bell’s histograms. The very good matchings at final flare 
and decrab were considered adequate validations of computations of distri- 
butions for these variables. 

C PTD 


With Bell’s control system, q)j D is not sensitive to the above mentioned 
drift situation that develops after decrabbing in a large crosswind. There- 
fore, we were able to obtain a meaningful comparison of bank angle distri- 
butions at touchdown. This comparison is shown in Fig. 28, where another 
excellent match is evident. 


TR- 1007-1 


70 



Figure 28. Comparison of Bell's Monte Carlo Generated Histogram of cprp-p 
■with the Computed Distribution (with Data for 
Both Guidance Systems Combined) 
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SECTION V 


CONCLUSIONS AND COMMENTS 


A. CONCLUSIONS 

The computed probability distributions were found to be extremely good 
fits to Bell's Monte Carlo generated histograms for pertinent longitudinal 
and lateral variables at two "points" along the vehicle's final approach 
trajectory. This validates the assumption that the longitudinal and lateral 
situations can be analyzed separately. As a further consequence of this 
good matching of results, it is apparent that (for the situation considered) 
the approximation technique has been validated as an alternative to a 
Monte Carlo simulation for obtaining distributions of pertinent variables 
through final approach and touchdown. However, due to the limited number 
of Monte Carlo runs available, this validation does not cover the extreme 
tails of the (computed) distributions. A significantly greater number of 
Monte Carlo runs would be required to validate computed points in the extreme 
tails of the distribution. However, as a practical matter, the analj^sis 
technique is most useful as a design tool in the region that has been vali- 
dated. Besides, the extreme tails of any "real" distribution are governed 
by factors that are not amenable to accurate representation anyway — either 
due to unknown inputs or to unknown extreme details of the distributions of 
known inputs. 

B. COMMENTS 

The principal advantage of the system model over a Monte Carlo alternative 
is its inherent flexibility. A Monte Carlo simulation will provide probability 
distributions, but system sensitivities are not direct outputs of such a 
simulation. If some aspect of the overall situation is modified, then a 
repeat of all the Monte Carlo runs is required to determine the effect of 
the modification. This is not only expensive, but it also fails to provide 
any insight into the fundamental interactions and tradeoffs that are present 
within any system. Such is not the case with the analysis technique and 
associated approach and landing system model. Touchdown sensitivities are 
an intermediate calculation to the estimated touchdown distributions. As a 
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consequence, the model, as used, is a highly practical analysis tool which 
obviates the need for a more expensive Monte Carlo simulation, previously 
required to obtain distributions or exceedance probabilities. 

It is pertinent to note here that several of the more detailed aspects 
of the analysis technique were not described in this report. Some of these 
aspects were omitted because they entail more or less standard mathematical 
procedures. Others were omitted because they were not used in the example 
applications presented. It is felt that some of these aspects deserve 
mention, if only to point out their existence. One area wherein several 
detailed facets of the analysis technique were omitted involves an appropriate 
representation of effects of random gusts after the start of flare. For 
example, it is possible to compute (and make use of) a time-varying rms value 
for dispersions during flare. Another point concerns the choice of flight 
path angle to be used in calculating the X T p sensitivity to h e at the start 
of flare. For low flare heights this isn't a critical matter. But for the 
higher flare heights (usually associated with one-step flares) the appropriate 
flight path angle appears to be that attained about 5 sec prior to touchdown, 
father than that during equilibrium glide . In addition, there are probably 
other detailed considerations where specific applications might require 
modification of the specific procedures presented here. 
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APPENDIX 


GUST INPUT TRANSFER FUNCTIONS JUST PRIOR 
TO FLARE (ONE- STEP FLARE)* 


.0206(. 0972)0 .lUO)(lU 73)[.847, 

. 0311 )(.1039)(.i493)(.- 458) (1.295] 


.206] [.71 5) 1-131 1[-7^1 > 4.32] 
;i4.io)L.523, - 928 J L-943, 3 . 9 V 


. OOl 48( . 993 ) ( 7 • 91 ) ( — 846 . ) [ . 999 » . 1O9U] [-.0986, .726]f.921, 1.398] 




